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DEVICE CONVERTING THERMAL ENERGY INTO KINETIC ONE 
BY USING SPONTANEOUS ISOTHERMAL GAS AGGREGATION 



My invention is a device converting thermal energy into kinetic one, related to 
the group of machines using four-phase basic thermodynamic processes like 
Carnot or Otto cycles. These devices need, for their operation, some kind of 
available outside heat source to be converted into kinetic energy. They consist 
of continuously lubricated moving parts, working in high temperatures, with 
quality deteriorating by usage and with noise emission. 
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My invention uses rarefied gas in a novel three-phase thermodynamic cycle, as 
shown in Fig.l (p,v) diagram, of which the first phase is a spontaneous 
isothermal gas aggregation (0 — 1), equivalent to an ideal isothermal 
compression, the second phase is an adiabatic expansion (1 — 2), with work 
produced via an expander and the third one is an isobaric expansion (2 — 0) 
where, by means of an exchanger, the cooled gas is reheated again (q2) by 
cooling the ambient air. The shaded area below the adiabatic path (1 — 2) 
represents the work done at the expense of the internal thermal energy of the 
gas (Iso). The first phase arises when the gas passes through numerous special 
microscopic holes, with sizes comparable to the mean free path of the 
molecules, so that the latter do not collide with each other but only with the 
walls. I have thought up smart geometric shapes for these holes, like slot 
(Fig.2) and cone (Fig.3) with diverging inner surfaces, cavity (Fig.4) with 
segments of spherical inner surfaces, in order that the molecules may take 
advantage of a phenomenon (to be discussed further down the text), with the 
result that, during successive rebounds upon the inner walls, they tend to move 
forward, forming a small but discrete net flow from the input (i) to the output 
(o). The solid lines with the arrows show the central paths of the swarms of 
molecules. Under these special conditions the gas comes out of the holes 
spontaneously and isothermally, entering a room with increased density. 
Obviously, there result five advantages, by the use of my invention, i.e. 
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(1) energy production at the expense of the internal thermal energy of the gas, 
which then is reheated by the ambient air, (2) refrigeration for any domestic 
appliances, (3) no moving parts (except the expander), (4) high quality 
operation and (5) no noise. 

5 

Description. 

Fig.5 (parallel view and cross section S-S) shows the device, consisting of a 
vacuum glass vessel (1) divided into two rooms (2) and (3) by a region (4) 
containing the microscopic holes assembly, consisting of a great number of 

10 holes, grouped into standard small modules (m), all arranged in a parallel 
layout as regards the gas flow. The closed circuit of the gas flow is 
supplemented with an adiabatic expander (5), within room (3) and a heat 
exchanger (6) in the return path of the gas from (3) to (2), transferring heat 
from the ambient air (7) to the gas with the help of ventilator (8). With suitable 

15 pressure difference between (2) and (3) an optimum flow is established^ so that 
the device is continuously performing work, e,g. by means of a generator (9), 
coupled to the expander through a magnetic clutch (10) and a speed reduction 
gear (1 1) (if needed) and at the same time it offers cooling possibilities. 

■ 

20 The Phenomenon. 

The operation of the device is based on a phenomenon observed at the time of 
the experimental research and evaluation of the external friction of gases [1], 
where it was shown that the molecules in a rarefied gaSj rebounded from the 
inner walls of the container, under suitable vacuum pressure, do not exactly 

25 obey the so called cosine-law (uniform rebound to all directions) [2,p.27], but, 
due to the existence of a molecular layer, adsorbed upon the walls, their path 
directions tend to slightly incline towards the perpendiculars to the walls, 
provided that the inner surfaces are quite smooth and the size of the container 
comparable with the mean free path of the molecules. Both of these properties 

30 are very important. The surface smoothness inside the holes must be perfect 
enough for the adsorption layer to cover the surface irregularities completely, 
otherwise the layer action is cancelled and the cosine-law prevails 
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again. Fortunately, nowadays a state-of-the-art value of surface roughness has 
been realized down to 1 nm, and even better [3], while in earlier decades 
values of less than 20 nm apparently had not been reached [4, p. 622]. With 
regard to the size, I have taken the fundamental dimension of the holes £ = 

^ 10 jam, which size is relatively easily realizable, happily in accordance with 
the technological progress of these days on Micro-Electro-Mechanical- 
Systems (MEMS) [5, p.56] and which is conveniently adaptable to the selected 
mean free path X = 10 jxm, as well as to the corresponding pressure [6, p. 24], 
within the range of a well developed molecular layer. Finally, I consider worth 

*® mentioning that this peculiar behavior of the molecular layers offers a natural 
explanation of the repulsive forces between adjacent corpuscles in the 
Brownian motion phenomenon and also in the expansion of dust in the air [1, 
p.331]. 

Industrial Applicability. 

The device has not been realized and tested experimentally. Nevertheless, its 
successful working ability is indeed proved indirectly, because it is based on 
the experimental and theoretical work mentioned in [1], as well as on a 
simulation method, assisted by electronic computer programs, to be described 
quantitavely as follows. 
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The Simulation Method. 

In order to evaluate the amount of flow through the microscopic holes, it is 
necessary first to calculate the number of molecules emitted from any point A 
of the inner walls and fallen on any other point B as a function of the 
geometric parameters (dimensions,- angles) of the holes. 
Following the computer symbolism, let 

AB[m] = distance between two points A and B located anywhere on the inner 

walls of a hole. 

*> 

na [sw/m ] = swarm of molecules per unit volume (volume density) around A 
dna [sw/(m *s)] = swarm of molecules per unit area per unit time rebounded 
from A within an infinitesimal stereo-angle dfi [sr] towards B. 
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v[m/s] = arithmetic mean velocity of the molecules 

cfa,cfb = cosines of angles <Pa,<Pb between AB and the perpendiculars on the 
respective infinitesimal facets dsa and dsb at A and B . 

na*v/4[sw/(m *s)] = molecules per unit area per unit time (surface density) 
5 rebounded from A to the inner hemisphere. 

Then, in the absence of the adsorbed layer the cosine-law is expressed as 
follows [2 ? p.27]j (Pi means 7t): 

dna = na*v/(4*Pi)*cfa*dQ = na*v/4*cfa*cfb/(Pi*AB 2 )*dsb 

Or, in reduced form (divided by no*v/4 and multiplied by dsa/dsb) 

10 dna*dsa/(no*v/4*dsb) = wa*cfa*cfb/(Pi*AB 2 )*dsa ( 1 ) 

where wa = (ria*v/4)/(rid*v/4) = relative surface density on A, wb=rio*v/4 = 
input surface density. On integration of dO over the inner hemisphere we 
obtain the basic quantity na*v/4. The factor cfa expresses the cosine-law. 
Now, in the presence of the adsorbed layer the cosine-law is to be modified, 

15 i.e. the factor cfa should be substituted by [1, p.325] 

{ [1-2/3 *f(p)]*cfa+f(p)*cfa 2 }, where f(p) is an increasing function with the 
pressure and with f(p) max = 3/2, occurring at p=l,918mmHg, which 
corresponds to (3/2* cfa 2 ) as a substitute of cfa. In that case 
dna*dsa/(no*v/4*dsb) = wa*3/2*cfa 2 *cfb/(Pi*AB 2 )*dsa ( 2 ) 

20 This formula may be used at least also for pressures above 1.918 [mmHg], up 
to 23,2 [mmHg], which corresponds to the maximum thickness of the layer 
and beyond, given that it does not drop quickly after the maximum 
[l,p.305,Table]. 

The forms of the holes are selected to possess some kind of symmetry so that 
25 the inner walls, as reflecting surfaces, may be divided into a large number (n) 
of strips (for the slots) and rings (for the cones and cavities), as shown in (12) 
of Figs 2,3,4. The same thing may be done on the input (i) and output (o) 
surfaces. Then, the relative density wa is constant along a strip or a ring. I 
have to remark that wa, when referred to the walls is an unknown, while when 
30 referred to the input surface it is known and equal to 1, and when referred to 
the output surface it is equal to the compression factor k between input and 
output. So, for each point B we are allowed to integrate (sum up) equations ( 1 ) 
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and (2) over each strip or ring, having previously expressed these equations as 
functions of the geometric parameters of the holes. After integration (addition) 
and by putting i for Ai(=u, 3 ,...n) and j for Bj(= u , 3 ,...n), I rewrite equations (1) and 
(2) in a new form 
swy = Wi*fbbpjj (layer absent) | 

swijij = Wi*fbbpij (layer present) | ( 3 ) 

where swy = swarm of molecules per strip or ring per unit time, rebounded 
from the strip or ring containing A\ to Bj per unit area for B. 
fbbpy = transmission coefficients from a strip or ring i to point j, that are 
calculated as functions of the geometric parameters. In order to find the n 
unknown densities, I express, in the form of equation, the following equality 
which, under steady-state conditions takes place between the number of 
molecules fallen on any reflecting point j and the number Wj rebounded from 
the same point. 

Si ( =u,3,...n)SWij[reflecting surface]+£ i( =i, 2,3,... n)SWy [input surface] 

+ k*Ei ( =i s 23,...n)Sw i j[output surface] wj ( 4 ) 

The first sum includes the unknown variables w*. The second and third sums 

are known. In terms of equations (3) this equality, appropriately rearranged, 

becomes an n-variable linear equation for point j : 

2i(=i,2,3, j-i) fbbp y *Wi + (fbbp ij -l)*w j + j+ 2 ,...n) fbbpij*Wi = 

- Zi ( =i,2 s 3,. n)fbbpij(input) - k*2i ( =i,2,3,..n)fbbpij(output) ( 5 ) 

Finally, we have a system of n n-variable linear equations, which may be 
solved with the help of Gauss algorithm [7,p.44-28]. 



Three examples. 

Having established the numerical values of the n variables (densities), both for 
layer absence and layer presence conditions, it is easy to calculate the 
algebraic sum Fl(k) of flow of molecules through the input or output (it is 
the same), including all the path combinations. This net overall flow Fl(k) is a 
linear function of k, reduced to the unit of input surface density no*v/4 and to 
the unit of area 1 Q 2 (slots and cones)[Figs 2,3] and r 2 (cavities)[Fig.4], (1 G = 
2*^, r = £). 
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Under layer absence and for k = 1 we have Fl(l) = 0* which complies with the 
cosine-law. Under layer presence and for k = 1 we have Fl(l) = Flm(maximum) 
and for k = km(maximum) the flow stops, i.e. Fl(km) = 0. Under layer presence: 
Fl(k) = Flm*(km-k)/(km- 1 ) ( 6 ) 

Flm and km are also functions of the geometric parameters of the holes; i.e. li;G> 
for slots and cones (Figs 2,3) and ac0,bd0 for cavities (Fig.4). Optimum values: 



Geometric parameters 
li(=li/lo) 
co[rad] 
ac0=bd0[rad] 
Overall flow Flm 



slot 
0,4 
1,4 



cone 
0,5 
0,8 



cavity 



0*052 



0 $ 0218 



0,7227 
0,1600 



Compression factor km 1,1100 1 ,2500 1 ,2000 

km is found by the trial-and-error method or directly with the formula: 
km = (A-Flm)/A (A=result of program with layer present,k=l,zero input) (7). 
Because of the great number of holes needed to achieve a somewhat remarkable 
result, I have organized the construction of the device in a form of small 
modules, as shown in Fig.6, consisting of a certain number (s) of parallel very 
thin panels, say xe(=0,3cm)*ye(=2,lcm), each perforated with a number of holes 
((13) for parallel slots of length all the way of the module's y-dimension, (14) 
for cones and cavities) and arranged in a pile(15) of height 
H(s) = s*h+2*d ( 8 ) 

where h(=0,2cm) = distance between successive panels, d(=lcm) = input or 
output air ducts. The arrows show the path of the molecules. Suitable supporting 
rods ((16), solid lines) fix the panels in place. Along z we have (s) holes in series 
and the molecule compression factor is k s (=kj *k2* . . . *k s ), (ki=k2=. . .=k s =k). 
The number Nmod(=ax*ay) of holes per panel or of piles of holes per module is 
estimated to Slot Cone Cavity 

Nmod = ax*ay = 80*(2cm/lo) 100*400 66*400 (9) 

Two gases, Helium and Hydrogen, have been chosen as the most suitable for use 
with the device. The present examples will work with Hydrogen (mass g[kg] = 
0,3347/1 0 26 , arithmetic mean velocity v[m/s] = 1693 [6,p.323] ). 
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Now, Fig. 7 (not in scale) shows a possible arrangement (18) of these modules 
(m) within a part O = 0,04241 m 3 (W=0,054 m) of a space (17) with dimensions 
X = lm and D (diameter) = lm, which will contain the device of Fig.5(modules' 
assembly and expander). I have taken a limited value of O in order to 
accommodate a heat exchanger of reasonable size for the device. The arrows 
indicate the gas flow directions (i=input, o=output). Then, the number v(s) of 
modules contained in O and the whole number Np(s) of piles of holes is, 
v(s) = 0/(xe*ye*H(s)) and Np(s) = Nmod*v(s) ( 10 ) 

With regard to Fig.l : Work done per cycle(shaded area) [8,p.244] 



(11) 



Ls [J/kg] = R [J/(kg*K)]*To[K]/ 

(n-l)*{l-(l/k s ) ((n - ,)/n) } 
R [4,p.872] = 4124, n [4,p.872] = 1,409 

To [K] = 253 for slots, 273 for cones and cavities (see next paragraph). 
In order to maximize the output power, the following expression a(k), which is a 
product of three factors in Eqs (6),(8),(ll),contained in the power output 
formula, must be maximized with respect to(k) and with (s) as a parameter, 
given that (s) may not exceed a limit (so),where the mean free path still remains 
"free" within the last holes, 

a(k)=(km-k)/(km-l)/(s*h+2*d)*{l-(l/k s ) ((nl)/n) } ( 12), 

to find k=ko, s=so. Computed values of ko, so, Fl(ko),H(so), v(so), Np(so), lso 
follow: slot cone cavity 

ko 1,05225 1,106 1,085 



so 

Fl(ko) 
H(so) [cm] 
v(so) 

Np(so)/10 6 
lso [J/kg] 



17 
0,0273 

5,4 
12465 
997,2 
566933 



0,01256 
3,8 
17715 
708,6 

637950 



11 

0,0920 

4,2 
16028 
423,1 

630466 



With plenty of margin (h) between successive panels and ample input-output ab- 
ducts (d), the speed of flow outside the holes is kept within a few meters per 
second, practically eliminating friction losses and noise. 
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Expander and Heat Exchanger . 

The expander [9, p.449] is a single - stage reaction gas turbine, accommodated 
within the device (Fig.5, (5)). Its main features of interest here are the wheel 
diameter (D), the revolving speed (n) and the efficiency factor pexp=0,825[9, 
05 p.271]. 

The exchanger[4,p.470-472] is constituted of 30 glass tubes (Fig.5,(6)) in 
parallel, 0,05m in diameter, lm of length, situated along and around the device; 
The gas H2 passes (in laminar flow) through the tubes, while £ir (Fig.5,(7)) is 
forced (in turbulent flow) around them, in the opposite direction, as shown by 

10 the arrows, by means of the ventilator (Fig.5,(8)), with velocities 2 to 5 m/s. 

In order to realize stich a reasonable size of this component, it was necessary to 
let a greater temperature drop between warm air and cool H2 (40°C for slots, 
20°C for cones and cavities). Fig. 8 shows schematically [9, p.271] the heat 
exchanger and the corresponding flow diagram. The horizontal and slanted 

* 5 arrows show air- and Hb-flow, vertical arrows show heat-flow. The (computed) 
pressure drop, in the Fb-flow is too small to be taken into consideration. 
Calculated values of (D), (n), and the working pressures and temperatures are as 
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follows (c v [kcal/(kg*K)] = 


2,41[4,p.871], e [kcal/J] = 0,2388/1 0 3 ) : 




Slot 


Cone 


Cavity 


EXPANDER D[m]-n[rev/min] 0,60-3630 


0,41-3630 


0,44-3630 


Pressure input pi=po*ko A so 


1020*2,377 


1121*2,48 


1121*2,45 


output po [Pa] 


1020 


1121 


1121 


Temperatue input To(=Td) 


253 


273 


273 


Output Tc=To-|3exp*lso*e/Cv 206,7 


220,8 


221,5 


EXCHANGER Input air tempTa 293 


293 


293 


Output air temp. Tb 246,7(-26,3°C) 


240,8(-32,2°C) 


241,5(-31,5°C) 


Input Hi temp. Tc 


206,7 


220,8 


221,5 


Output H2 temp. Td(=To) 


253 


273 


273 


Ta-Tb =Td-Tc 


46,3 


52,2 


51,5 


Air flow rate[m /s] 


0,95 


0,66 


0,77 


Ventilator Power I vent, [w] 


190 


120 


140 
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Hydrogen reheating thermal energy (Fig.l) [8, p:2-35]:q2= e p *(To-Te) 
(c p [kcal/(kg*K)] = 3,40[4,p.871] Slot Cone Cavity 

q 2 [kcal/kg] 157,42 177,48 175.10 

Numerical Results. 

Finally; I proceed to calculate all the factors which determine the output power: 
Loschmidt number [6, p. 17] (p=l,02*10 5 Pa, T=273K ) = 

= 2,687* 1 0 25 molecules/m 3 
Input pressure po[Pa] 

po[mmHg] 
Input Temperature To[K] 
Input Vol. Density no [sw/m 3 ]/10 23 
Hydrogen Velocity v [m/s] 
Input Surf. Density: 
wo = (no*v/4) [sw/(m 2 *s]/10 23 
lo [m] = 20/1 0 6 r[m] = 1 0/1 0 6 

Mass flow rate per hole: Slots and Cones 

Cavities 

Total flow rate 
Power output of expander 
Power output (pract.) 



Slot 


Cone 


Cavity 


1020 


1121 


1121 


7,68 


8,41 


8,41 


253 


273 


273 


2,900 


2,950 


2 5 950 


1630 


1693 


1693 



1182 



1249 



1249 



gf[kg/s] = g*Fl(ko)*wo*lo 2 
gf[kg/s] = g*Fl(ko)*wo*r 
G[kg/s] = gf*Np(so) 
Iexp[watt] = ftexp*lso*G: 
Ipr[watt] = Iexp - Ivent 





Slots 


Cones 


Cavities 


Fl(ko) 


0,0273 


0,01256 


0,0920 


g frkg/s]*io 12 


4,32 


2; 10 


3;85 


G[kg/s]*10 3 


4,308 


1,487 


1,629 


lsb[J/kg] 


566933 


637950 


630466 


Iexp [watt] 


2015 


783 


849 


Ivent[watfJ 


190 


120 


140 


Ipract[watt] 


1825 


663 


709 



Construction Hints. 

Mass production can be achieved by the method of pressing [10,p. 8-1], 
not excluding any other competent method. As construction material I would 
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propose glass, ceramic, silicon or the like, used in semiconductor technology. 
Fig.9 shows a slot panel i.e. an arrangement of parallel triangular rods (19), 
forming slots (s) in between, lying on supporting rods (20) (cross-section Ti-Ti) 
at suitable intervals. Cross section T2-T2 of rods (19). The distance between 
successive panels is h=0;2 [cm]; Both forms of rods can easily be manufactured 
in mass production with the active surface (b) made very smooth by advanced 
polishing processes [5* p i 56] i 

The slot solution presents evident advantages over the other two solutions in, (a) 
manufacture, (b) greater output power per unit volume. 

Fig. 10 shows a cone panel (21) with cones (e) (cross section T2-T2), arranged in 
series along x, lying on supporting rods (22) (cross-section T1-T1), which are 
placed between adjacent cone series. Intervals between successive panels are 
equal to h=0,2 [cm] . The cone active surface (b) is made very smooth. Fig. 1 1 
shows a possible scheme for cone panel fabrication* with the help of molds 
(2a,cylinders), (2b) and (p) as pressing means. 

Finally, Fig. 12 shows a cavity panel (23), carrying the holes with the active 
spherical surfaces (b) and the supporting rods (24) (cross sections (Ti - Ti , 
T2- T2 ) ), carrying the active spherical surfaces (c). At suitable intervals along 
the rbds(24), a contact rod (25) is made in place of the corresponding active 
surface (c), with elimination of the opposite side hole, in order that the panel is 
rigidly supported. Figs 1 3 and 1 4 show the forming of the active surfaces (b) and 
(c) of the cavity respectively, with the help of molds (3a),(3b),(3c,cylinders),(p) 
for Fig.T3 and (4a),(4b)*(p) for Fig; 14; To achieve the exact spherical surface the 
molds should be equipped with tiny balls (s, dia. 20jim), with smooth spherical 
shape, ike those used in miniature ball-bearings [11]. 
Computer Programs. 

A 3 V2 in. floppy disc is available^ containing the programs (written in Q-basic) 
of the present invention. 
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Brief description of drawings (amended) 



Page 1/7 (amended). 

Fig. 1 : Three - phase thermodynamic cycle diagram. 

Fig. 2: Slot in face and side view. Diverging inner surfaces. Microscopic cross section 
(side view). Macroscopic length 20mm (30). 
Fig. 3 : Cone in face and side view. 
Fig. 4: Cavity in face and side view. 

Page 2/7. 

Fig. 5: Complete device of the invention (side view and cross section S-S). 
Page 3/7 (amended). 

Fig. 6: Small modules (m), panel (13) of parallel slots, cone and cavity panels (14), 
pile of panels (15). 

Fig. 7: Arrangement (18) of modules (m). The arrows indicate the directions of the 
gas stream. 

Page 4/7. 

Fig. 8: Heat exchanger diagram, corresponding flow diagram. 
Page 5/7. 

Fig. 9: Slot panel (arrangement of parallel triangular rods (19), forming slots (s) in 
between). 

Fig. 10: Cone panel (21). 
Page 6/7. 

Fig.l 1 : Possible scheme for cone panel fabrication. 
Fig. 12: Cavity panel (23). 
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Fig. 13: Forming of active surfaces on cavities. 
Fig. 14: Ditto 



